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Nonlinear hydrodynamic equations for visco-elastic media are discussed. We start from the recently 
derived fully hydrodynamic nonlinear description of permanent elasticity that utilizes the (Eulerian) 
strain tensor. The reversible quadratic nonlinearities in the strain tensor dynamics are of the 'lower 
convected' type, unambiguously. Replacing the (often neglected) strain diffusion by a relaxation 
of the strain as a minimal ingredient, a generalized hydrodynamic description of viscoelasticity is 
obtained. This can be used to get a nonlinear dynamic equation for the stress tensor (sometimes 
called constitutive equation) in terms of a power series in the variables. The form of this equation 
and in particular the form of the nonlinear convective term is not universal but depends on various 
material parameters. A comparison with existing phenomenological models is given. In particular 
we discuss how these ad-hoc models fit into the hydrodynamic description and where the various 
non-Newtonian contributions are coming from. 

PACS numbers: 05.70.Ln, 46.05.+b, 83.10.Nn 

I. INTRODUCTION 

Hydrodynamics is a well established field to describe macroscopically simple fluids by means of the Navier-Stokes-, 
continuity, and heat conduction equations. However, it applies also to more complex fluids that arc fully characterized 
by conservation laws and broken symmetries. This more general hydrodynamic method has been established in the 
60s [B 12- 13 and applied e.g. to superfluids and liquid crystals [EJ. It is based on (the Gibbsian formulation of) 
thermodynamics |6J, [jj , symmetries and well-founded physical principles 8J . A detailed description of this method can 
be found in 0, . Somewhat related approaches have been used for liquid crystals [ToL ITTl \VA , IT^ | and more generally 
in [0 C3 

On the other hand, non-Newtonian fluids are believed to show non-universal behavior and a host of different em- 
pirical models have been proposed 0, 0, 0, 0, |2(1 to cope with the flow rheology of such substances. Typically 
these models are formulated as generalizations of the linear, Newtonian relation between stress and deformational 
flow allowing for additional time derivatives and nonlinearities. They are tailored to accommodate empiric findings 
or are based on principles |20| that are ad-hoc and generally insufficient. 

Quite recently we have derived a nonlinear hydrodynamic description of elastic media |22l |23| that is based on 
first principles, only, making use of thermostatics, linear irreversible thermodynamics, symmetries and broken sym- 
metries, and invariance principles. It has been confirmed within the GENERIC formalism Allowing in this 
hydrodynamic description the strains to relax (and not only to diffuse) a generalized hydrodynamic description of 
nonlinear viscoelasticity is obtained in terms of a dynamic equation for the strain tensor p2l |23| . After a thorough 
exposition of this formulation (Sec[Hj we transform it into a description in terms of a dynamic equation for the stress 
tensor (Sec lIII|l . This can only be done approximately in the form of a power expansion in the variables. Up to 
second order a formulation is obtained that can directly be compared with many of the empirical models proposed to 
describe non-Newtonian rheology. The comparison fSec l! V[) reveals possible inconsistencies and connects the various 
ad-hoc additions of those models with physical relevant processes, like strain relaxation, elasticity and viscosity. In 
two appendices we sketch possible extensions of the hydrodynamic description of viscoelasticity. 

A comparison with recent constitutive equations that refer to specific microscopic variables and processes, like 
convective constraint release [2^, [2(| [27j , will not be done in the present paper. Here we rather concentrate on the 
simplest generalized hydrodynamic description of non-Newtonian rheology in terms of a relaxing strain field, while a 
detailed comparison with those theories requires the use of additional relaxing fields. 
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II. STRAIN TENSOR DESCRIPTION 



In this section we review the hydrodynamic description of nonlinear elasticity and its generalization to viscoelasticity 
in terms of the Eulerian strain tensor Uij. We start with a short reminder of its definition. In a stress-free elastic 
body, we consider a point at the (initial) coordinate a. As the body is displaced, rotated, compressed and sheared, the 
given point is displaced to r. Since all points of the body have a unique pair of a and r, the function r(a) is unique and 
invertible, the result of which we denote as a(r). Describing physics in local terms, in particular, the state variables 
characterizing the elastic body are taken as functions of the real space coordinates r rather than of the reference space 
coordinates a. Therefore we choose a(r) rather than r(a) as the starting point to construct the (Eulerian) strain tensor 
Uij — h[8ij — (Vj£ife)(Vjafe)], which can be written in the more familiar form Uij — ^[Vj-u; + V^Mj — (Vjitfe)(Vjitfe)] 
using the displacement field Ui(r) = ri — di(r). In the Lagrangian description, where all fields are function of the 
initial coordinates rather than the local ones, the (Langrangian) strain tensor is U^j — \[{drk/ ' da,j)(drk / dai) — <%] 
or Uij = ^[dui/daj + duj/dai + (duk/dcii) ■ (duk/daj)], with Ui(a) = ri(a) — a,. Since the dynamics for a is simply 
da/dt = in the absence of phenomenological currents, the dynamics for the Eulerian strain tensor reads p2l |23| 

d 
dt 



^Uij - Aij + U ki VjV k + U^m = x\f ] (1) 



with d/dt = d/dt + v • V, where Vi is the velocity and = i(Vji>j + Vjt>i) its symmetrized gradient describing 
deformational flow. This equation contains, apart from the transport derivative ~ v • V due to Galilean translational 
invariance, a linear and a nonlinear coupling to flow. The former (~ A t j) reflects the spontaneously broken transla- 
tional symmetry in elastic media (in an equation for the displacement vector the equivalent term subtracts the solid 
body translations, which must not change the elastic state of the fluid). The nonlinear one is of the lower convected 
type and results from the freedom to choose the orientation of a material-fixed frame independently from that of the 
laboratory frame |22l \2?fy . For a Lagrangian strain tensor the upper convected derivative would be obtained. On 
the right hand side there is the phenomenological (quasi-) current X^ h \ which describes, in the case of permanent 
elasticity, purely diffusional processes. In that case Eq.Q), when combined with the dynamic equations for the other 
hydrodynamic variables, the density, p, the energy density, e, and the density of linear momentum, gi, gives rise 
to three additional truly hydrodynamic modes, i.e. the two doubly degenerated transverse sound modes (compared 
to the two vorticity diffusions of simple fluids) and vacancy diffusion |5| . For decaying strains in viscoelastic media 

X^j h ^ contains relaxational processes (see below). The strain tensor Uij is symmetric, in order to exclude solid body 
rotations, since the latter must not change the elastic state. 

The second relevant dynamic equation is the momentum balance or generalized Navier-Stokes equation 

d 

p-gui + V lP + Vj<j i3 = (2) 



with the isotropic pressure p and the stress tensor , which is of the form 



<Tij = -*« + VuUjk + * fe j Uik + (3) 

Here ^ij is the elastic stress tensor, the thermodynamic conjugate to the strain tensor. The linear and nonlinear 'Fy 
contributions are the counterterms to the linear and nonlinear flow contributions in Eq. Q that are necessary to 
cancel any contribution to the entropy production |9|, since all these terms are reversible. Thus, these parts of the 
stress tensor are completely fixed by general physical principles. The phenomenological part of the stress tensor is 
contained in <J~^ (see below) and describes in the simplest form Newtonian viscosity. The stress tensor is symmetric 
(or can be made so) in order to guarantee angular momentum conservation [{|. 

Throughout this paper we will assume incompressibility that reduces the continuity or mass conservation equation 
to divu = 0. In addition we neglect the thermal degree of freedom and, thus, the heat conduction equation. 

Phenomenological or material specific properties enter our equations on three different occasions, two are dynamic 
and one is static in nature. For the former we stay within the reach of well-founded linear irreversible thermodynamics 
psj . i.e. we assume a linear relation between the "forces" ^fc; and Aki and the "fluxes" X^^ and cf\^ 

X^ = -a m ^u (4) 
= -VijkiAu (5) 

describing strain relaxation and viscosity, respectively. By allowing nonlinear dependences of the 'fluxes" on the 
" forces" , one would leave the solid grounds of well-established statistical physics, since not very much is known on 
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the validity range of such theories. Nevertheless, the relations (0] |3J are nonlinear in the sense that the material 
tensors depend on the variables of the systems, in particular on temperature T, pressure and the strain tensor. The 
dependence on the scalar quantities is rather trivial and will not be shown explicitly, while the strain tensor dependence 
is more complicated. Starting from an equilibrium state = we can expand the material tensors into powers of 
Uij. Up to quadratic order in Eas. (|1I2|I . which is what we need for the comparison in Chapt. IV, we find for the 
general form of the rank 4 material tensors 

CL\ CL 2 

a-ijki = -^(dikSji + SuSjk) + -^(UikSji + UjkSu + UuSjk + UjiSik) + a 3 SijS k i + 

ai{5ijUki + SkiUij) + a 5 S i: jSkiUpp + a e (Si k Sji + 8 a 5jk)U pp . (6) 

Since the " fluxes" are derived from a potential (the entropy production) by partial derivation, the symmetries of the 
tensors A+j and are transferred to <J^ and X^ h \ respectively and are reflected in the material tensors being 
symmetric under the interchange of i — j, k — I, and ij — kl. For the viscosity tensor Vijki we can simplify the form 
lO, since Akk = 0. Thus, all parts ~ 6ij or Ski are projected out in the entropy production and only the coefficients 
^1,2,6 appear in Eq.©. Note that neither Uij, nor ^ij, nor the total stress tensor er^ have to be traceless, despite Ay- 
being so. If desired, the expansion © can easily be continued to arbitrary order. A more general form of Eas. (|4l5|l is 
discussed in Appendix A. In the main part of this paper we will stick to the simplest form given above that constitutes 
the minimal generalization of hydrodynamics, but still captures viscoelasticity. 

The static phenomenological properties are expressed by a relation between the strain and the elastic stress 

^ij=K ijH U H (7) 

where the elastic tensor Kij k i generally depends on all variables, in particular on the strain tensor itself. In second 
order its form is slightly different from that of the a tensor, since it is derived from an elastic energy, which is cubic 
in the strain tensor, resulting in the constraint 2K§ = K4 

K^ki = -^-(SikSji + SuSjk) + -^-(UikSji + UjkSu + UuSjk + UjiSik) + K 3 6ijS k i + 



K^SijUki + SkiU^ + -[SikSji + SuSjk}Upp) + K 5 S t jSkiU P p. (8) 



Thus we are left with two linear elastic moduli K\ 3 and three quadratic ones ^2,4.5- To ensure positivity of the 
elastic energy, from which J7J is derived, third order terms with large and positive elastic moduli are necessary (see 
Appendix A). 



III. STRESS TENSOR DESCRIPTION 



In the preceding section viscoelastic hydrodynamics was expressed by the strain tensor and its derivatives. We will 
now rewrite these equations with the stress tensor as variable by replacing the strain tensor (and its derivatives) by 
the stress tensor (and its derivatives). This can only be done in an approximate way, since the equations are nonlinear. 
We will set up a power series expansion up to second order in the (old and new) variables. Of course, the resulting 
equations are less general than the starting ones and only applicable, if quadratic nonlinearities are sufficient for the 
problem at hand. It is then also sufficient to restrict the phenomenological expansions (|4I5I7|I to quadratic order. 
In addition, we will neglect in this section also those phenomenological terms that are connected with the traces of 
the variables (Ukk, ^fefc, ^kk )■> ^- e - there are only six static and dynamic material parameters left (ai,2, ^1.2, and 
.Ki^). A more general treatment including the trace-related second order phenomenological constants will be given 
in Appendix B. With this proviso the dynamic strain equation reads 

^-U i3 - A^ + U ki VjV k + Uk 3 ViV k = -—U^ - —U ik U 3 k (9) 
at t\ t 2 

with the abbreviations T-j - = a\K\ and t^ 1 = 2aiK 2 + 2a 2 Ki. The stress tensor that enters the Navier-Stokes 
equation @ is written as 

a U = ~ K i u ij + K 2 U ikUjk - viAij ~ v 2 {U ik A jk + U jk A ik ) (10) 

with K' 2 = 2K\ — 2K 2 . Linearizing in the phenomenological expansions (|4I5I7|I . i.e. K 2 = v 2 = t 2 x = 0, does not 
linearize the stress tensor, because of the convective nonlinearities in the strain relaxation equation 0. Only the 



4 



condition K 2 = (and v 2 — 0) linearizes the stress tensor, which is however a rather artificial condition, since it 
requires the linear and quadratic elastic moduli to be equal exactly. 

Taking the derivative d/dt of in l|10[l and replacing dUij/dt according to Eq.© we get 



j t <rij=-f[U ij ,A ij ,j t A ij ,n ij ] (11) 



in terms of the strain tensor, the deformational flow and its derivative, and the vorticity fiy = ^(\7jVi — V^u,) 
describing rotational flow. The latter enters due to the convective nonlinearities in 0. To convert this into the 
desired dynamic equation for the stress tensor we have to invert = aij[Uij, Ay], Ea. (|10() . into Uij = Uij[<Jij, Ay]. 
This is done approximately by the power expansion J7y — uf^ + [/(? uod ) + . . ., where U-j and u^ uad ^> contain 
expressions linear and quadratic in the variables, respectively. In particular we find 

KxU^ n) = (12) 
KfU^ uad) = K' 2 a lk a ok + (K' 2 v x + K lV2 ){a lk A ]k + a jk A ik ) + v x {K' 2 v x + 2K x v 2 )A lk A, ]k (13) 

Using these expressions the dynamic equation for the stress tensor takes the final form 

D a Df, r 

Ti-prr&ij + 0y = -Uoo Aij - ViTi— Aij + —o ik o ]k 

L)t L)t A] 

t v ( d \ 

+ V' 7jk + UlA ^di Aik + + U ^ A ik]g^ A jk) + 0(3) (14) 

where =V\ + T\K\ and r = ti/t 2 — K' 2 jK\ and 

= ~ a{Ti k Aj k + Tj k A ik ) — (T ik Qj k + Tj k Vl ik ) (15) 

for any tensor and number a. For a = —1 (a = +1) D a /Dt is the lower (upper) convected derivative, for a = 
the Jaumann or corotational derivative, while for a general a a linear combination of those is invoked. In our case the 
numbers a and b are 



v\ K' 2 v a 



K x t 2 K{ n 

v\ K 2 Vqq a 2 v 2 



(16) 



b = -l + — (17) 



In the very special 'linear' case discussed above (K 2 = t 2 x = v 2 = 0) we get a = b = — 1 and both time derivatives 
in Eq. (|14f> are of the lower convected type (as was the case for Uij without any approximation). However, this case is 
rather artificial and, generally, the phenomenological terms push these time derivatives away from the lower convected 
type and make them material dependent. 

Splitting Ea. (|14|l into the trace and the traceless part, it is easy to see that a kk and (d/dt)<r kk are of 0(2), since 
A kk = 0, and do not influence the dynamics of the traceless part cr? up to second order. This statement remains 
true, even if the trace-related second order phenomenological constants are taken into account (Appendix B). 



IV. COMPARISON WITH CONSTITUTIVE MODELS 



Ea. (|14|l is in a form that many empirical constitutive models have and a direct comparison is possible. First there 
are the models of the form 

nj^Vij+Vij = -VooAij (18) 

that lack the time derivative of Ay and the terms ~ r and ~ v 2 . For a — —1,0, 1 these are the Maxwell models 
poj . while for —1 < a < 1 it is the Johnson-Segalman model |29j . Apparently in these models the Newtonian 
viscosity v± has been neglected. For the actual viscosity that means it is given by t\K\ alone. This seems to be 
a good approximation for polymer solutions, where v\ is interpreted as the viscosity of the solvent and t\K\ as due 
to the polymers. Of course, if v\ = 0, it only makes sense to also put v 2 = 0. In addition, r — generally implies 
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t 2 x = = K' 2 (the relation ti/t 2 = K' 2 / Ki that also leads to r = is rather singular) meaning that these models are 
all within the linear case described above. However, this in turn has the consequence that a = — 1, which makes the 
Johnson-Segalman model inconsistent. The Maxwell model with a = 1 is applicable in a Lagrange description, while 
a = is not appropriate for a stress tensor like quantity (but rather for an orientational order parameter .30]). 
The next group of models take into account the time derivative of Ay 



Da A , D a 

— — c 
Dt 



r i7^ ff y + a >i = —Voo-Aij - z^oo-V-p^Ay (19) 



but neglect the terms ~ r and ~ v 2 . Thus, they again are in the linear class, a fact which is compatible with a = b 
that is inherent to these models. Again a = — 1 (Oldroyd A) is required for the Eulerian case (and a = 1, or Oldroyd 
B, for the Lagrangian description) while a = (Jeffries) is inappropriate. 
The first model that goes beyond the linear case is the Giesekus model 

Di 2cti 

Tl—(Tij+(Tij = -VooAij OikOjk (20) 

Dt Voe 

with the phcnomenological number c with < c < 1. This model neglects v\ and consistently also v 2 . Since c is 
positive, comparison with l|14(l shows that r < or K 2 jK\ > t\/t% is assumed. However, with K' 2 ^ one immediately 
concludes from Ijltjfl that a has to deviate from the pure value —1 (or 1). Thus, a model with nonlinearities in the 
stress tensor like the Giesekus model must not be of the pure lower (or upper) convected type, but has to show a 
material dependent deviation from that. For v\ = = v 2 this deviation from a = — 1 is — 2c/ Ki. 

All the models above are restricted to nonlinearities quadratic in the variables and do not show any term we have 
hidden in 0(3) in Ea. (|14(l . This is different in the PTT and PT models j^lj,l22j that have a stress tensor nonlinearity 
of the form caijUkk + 0(4), which for A k k — is of third order. These models have consistently a stress convection 
with a different from — 1. However, there are several other 3rd order terms (cf. Ea. ljB.7| and even 2nd order ones, 
not considered in these models. 

Another popular rheological model is the 'second order fluid' [2(j that contains a contribution 

o-.jf = -ci^^-Ay - c 2 A ik A jk (21) 

in the stress tensor (± for the Lagrangian and Eulerian picture, respectively). Comparison with Ea. (|10fl reveals that 
it is contained there with c\ = v\t\ and C2 = fiTi(6 + l) (for the Eulerian case), where b is defined in Ea. l|17|) . In the 
linear case defined above b = — 1 and therefore C2 = 0. Thus, the quadratic flow contributions to the stress tensor are 
due to the strain dependence of elasticity, viscosity and strain relaxation, when the strain field is taken into account 
as a relaxing variable and then eliminated. 

Even when no convective derivatives arc considered, very often nonlinear phcnomenological stress/strain relations 
qip h ) — F(A) are used. As an example, in a power law fluid the Newtonian shear viscosity v\ is replaced by 
v = Yl l n=i v n{^ikAik)^ n ~ 1 ^ 2 ■ Much more complicated forms are used [3j|. The problem with all these models is 
the compliance with thermodynamics. The expansion © for the viscosity tensor in terms of fTy avoids this problem 
and can be carried on to any order desired. However, similar contributions to the stress tensor originate from the 
expansions of a.ijki and Kijki- In addition, these expansions also change quite considerably the structure of the 
dynamic equation for the stress tensor (|14f> rendering inconsistent any model that uses a power law description of the 
shear viscosity, only. 



V. SUMMARY 



In this manuscript we have shown that the hydrodynamically derived model for non-Newtonian fluids in terms of 
the Eulerian strain tensor contains most of the standard rheological models as special cases and discards a few of them. 
Of course, the former is more general, as it contains powers of the relevant fields of arbitrary order when written in 
terms of the stress tensor. The hydrodynamic method allows to discriminate those parts of the dynamics that are due 
to general principles from the unavoidable phenomenological part. The latter is given here in the form of truncated 
power series in the strain tensor that can systematically be generalized when necessary. For the phenomenological 
part we stick to the well-established 'linear irreversible thermodynamics', which, being linear in the generalized forces, 
nevertheless leads to equations highly nonlinear in the variables like the strain tensor. 
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APPENDIX A: Generalization of the phenomenological equations 

In Eas. i|4l5|) we have omitted possible reversible crosscouplings between flow and strain dynamics 

X\f ] = -a ijk iV kl -p ijk iA kl (A.l) 
4j H) = -^jkiAki+PkUj^M (A.2) 

characterized by the tensor Pijki, which is of a slightly more complicated form than a^-jy in Eq.©. Being reversible 
it is not derived from a potential and lacks the ij — kl symmetry. This results in two different components (3<ia&ijU k i + 
Pib^kiUij^ where, however, the latter (as well as f3 3 ,s) vanish due to incompressibility. Thus we are left with 4 
parameters /?i,2.4a,6- Such terms, coming with probably small reactive transport parameters, have to compete with 
the parameter free, symmetry required terms already being part of Eas. i|ll2l) . These reversible crosscouplings are 
possible in the viscoelastic case only, and are absent for permanent elasticity. In the latter case only diffusion 
x (ph) = oVfefVi^ + Vj^ ik ) is present. 

We also list here the 5 third order static elastic contributions that follow from an elastic energy quartic in the strain 
tensor 

K ljki = K 6 Sij5kiU pp U qq + K 7 ([S lk Sji + Sudjk]UppU qq + dijSkiUpqUpq) (A. 3) 

+K 8 (SijU k pUip + 5 k iU ip Ujp + j[S l kU J i + S jk U l i + SuUjk + SjiUikjUpp) 
+K 9 (Si k Sji + 5u5j k )Up q U pq + Kio(8i k UjpUi p + Sj k Ui P Ui p + 5uUj P U kp + SjiUi P U kp ) 

and that are a generalization of the first and second order terms kept in Sec^ A necessary stability condition 
involving the second order modulus K 2 requires 27Kx(Kg + Kxq) > 2K 2 . 



APPENDIX B: Influence of the tensor traces 



Here we discuss the influence of those phenomenological parameters in Eas. <|6l7|l that are connected with the 
traces of the tensors involved, (23,4,5,61 v§, and -^3,4,5, as well as /3i,2.4a,6 <|A.1IA.2|) neglected in the main text. The 
phenomenological parts of the dynamic and static equations then read 

X^ h) = -a^j - a 2 (U ik ^j k + U jk ^ ik ) - a 3 S %3 ^ kk - a 4 {5 l3 U k i^ k i + U^ kk ) - a 5 5 l3 U kk ^ u 

-2a e U kk y i:j - pxA^ - /3 2 (U ik A jk + U jk A ik ) - f3 ia d t3 U k iA k i - 2(3 6 U kk A 2j (B.l) 

a ^ h) = ~^ A ij - "2(U lk A jk + U jk A lk ) ~ 2v & A ij U kk + px^ij + foi^ikUjk + ^ 3 kU tk ) 

+p4aVkkUij + 2As%f7« (B.2) 
*y = KxUij + 2K 2 U ik U jk + K 3 S l3 U kk + K A (5 tJ U kl U kl + 2U l3 U kk ) + K 5 S ZJ U kk U u (B.3) 

The strain relaxation and the stress take the form 

-jjrUij - Aij = -r^Uij - T 2 1 U lk V 3k - T 3 ~ 1 5 lj U kk - T^ l 8ijUkiUki - T 5 ~ 1 8 ij U kk Uu (B.4) 

— t 6 1 UijU kk — PxAij — f32{Ui k Aj k + Uj k A ik ) — P4 a SijU k iA k i — 2(3 e lI kk Aij 
= -KxUij + K' 2 U lk U jk - K 3 8ijU kk + K 3 UijU kk - K A 5 tj U k iU k i (B.5) 
-K' 5 5ijU kk Un - uxAij - v 2 (U ik Aj k + U jk A ik ) - 2v & A i] U kk 

with the new parameters T3" 1 = axK 3 + a 3 K\ + 3a 3 K 3 , t^ 1 = oi\K± + 2a 3 K 2 + 303^4 + a^Kx, t^ 1 = 2a 3 K 4 + 
a i K 3 + a 5 K 1 +a 1 K 5 + 3a 3 K 5 + 3a 5 K 3 + 2a 6 K 3 , r 6 _1 = 2a 2 K 3 + a i K 1 + 2a 1 K i + 2,a i K 3 + 2a & K u K x = K^l-fix), 
K 3 = K 3 {1 - A), K 4 = K 4 (l - fa), K' 2 =K' 2 + 2p 2 K, + 2PxK 2 , K' 3 = 2K 3 (1 + f3 2 + §/3 4a ) - 2K A + Kx(/3 ia + 2fo), 
andK' 5 =K 5 {l-Px)-2K 3 p 6 . 
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The procedure to switch from the strain relaxation equation to an effective dynamic equation for the stress tensor is 
the same as described in Sec. IIIII Due to the many new terms, it is more involved, and it is complicated additionally 
by the fact that the traces and the traceless parts of Uij and Uij are related by different parameters (even in the linear 
case) as can be seen from the modified Eg. <fP2|) 

KxUf^ = -<7 tj - v-iAij H 3 - 8ijcr k k (B.6) 

3^3 + K\ 

We refrain from writing down more details, but discuss the structure of the final result 
D a D~ f d 

t is ( d d 

+ -g— yP'jk + B 3 6 jk crpp + viA jk \—A ik + [a. lk + B^S^aqq + viA ik ]—A jk 

+ a- kk (B A a i: j + B 5 Atj) + S i: j(B 6 a k ia k i + B 7 a H A k i + B 8 A kl A kl + B 9 a kk a u ) + 0(3) 

with z/qo = v\ + KiT\(l — fii) 2 and the tilted numbers a, b, and f being much more complicated than the untilted 
ones I|16ll7fl . There are structurally new terms related to a kk , a linear one (with B\ = K 3 /(3K3 + K\) + t\/t%), 
and eight quadratic ones characterized by coefficients -B2,...,9- Nevertheless, the trace of the stress tensor and its time 
derivative are of second order and do not influence the dynamics of the traceless part in 0(2). Thus, for the deviator, 
afj — Uij — (l/3)<5y(Tfefc, the dynamic equation has exactly the form l)14|l . but with the tilted numbers and parameters 
instead of the untilted ones. For the trace we find to second order 

Ti4°fcfc + ( 1 + 3B iKfc = ^v 1 T 1 A k iA kl + ^al l al l + 2 ^^{al l +v 1 A k i)^-A kl 
at Ki Ki at 

+ 3(B 6 a Q kl a° kl + [B 7 + ^ Tl a] a° kl A kl + B s A kl A kl ) + 0(3) (B.8) 



indicating that its relaxational dynamics is completely given by the flow Aij and the traceless part a. 
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